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UNVEILING THE PHENOMENON OF DOUBLE PERIODIC VARIABLES 
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R.E. Mennickent 1 and Z. Kolaczkowski 1,2 
RESUMEN 

En este artfculo damos un reporte breve de nuestra investigation sobre Variables de Periodo Doble (DPVs), 
incluyendo el descubrimiento de DPVs en nuestra galaxia y algunas ideas acerca de la naturaleza de su periodo 
largo. 

ABSTRACT 

In this paper we give a brief report of our recent research on Double Periodic Variables (DPVs) , including the 
discovery of DPVs in the Galaxy and some insights on the nature of their long-cycle variability. 
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Double Periodic Variables (DPVs) are blue stars 
characterized by a short periodicity (1-16 days) and 
a long periodicity (50-600 days) in their light curves. 
They were discovered in the Magellanic Clouds after 
a search for Be stars in the OGLE variable star cata- 
log (Mennickent et al. 2003) and now they amount to 
« 110 in the MCs and about a dozen in our Galaxy 
(Fig. 1). A key feature of DPVs is the relation 
between short and long period, hereafter called su- 
perperiod. While the short period reflects the or- 
bital motion of a binary (Mennickent et al. 2005), 
the superperiod was an enigma until recently, when 
evidence was presented for an interpretation based 
on cycles of formation and dissipation of circumbi- 
nary discs (Mennickent et al. 2008; M08). We have 
found a significant infrared excess for all DPVs with 
available 2MASS data. To explain the DPV phe- 
nomenon, we have developed an hypothesis. In our 
view, DPVs are low mass ratio, semi-detached, in- 
termediate mass binaries in a stage of mass trans- 
fer, probably at the end of Case-A or Case-B mass 
exchange (M08). The B-type primary is rotating 
very fast while accretes matter from the later type 
donor in an Algol-like configuration. Matter is also 
ejected from the binary system through the outer 
Lagrangian points forming a reservoir of circumbi- 
nary gas. The superperiod could be produced by a 
variable mass ejection rate linked to the binary pe- 
riod. We have initiated an observing campaign to 
determine fundamental parameters of DPVs, mainly 
based on the study of the bright eclipsing systems, in 
order to shed light on the DPV evolution stage and 
their connection with classical Algol stars. In this 
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program, the sample of Galactic DPVs HD 50526, 
AU Mon, V 393 Sco and V356 Sgr are being stud- 
ied intensively, in optical and infrared wavelengths 
(Fig. 2). We have identified, based on high- resolution 
spectroscopy of some of these systems, at least three 
spectral components: (1) the secondary star revealed 
in absorption lines characteristics of a B-A type spec- 
trum, (2) an optically thin Ha emitting region (we 
observe sometimes double peak emission and some- 
times a single emission at one of the wings of the 
Balmer absorptions), and (3) an optically thick re- 
gion forming variable and broad helium absorption 
lines that could be signatures of the inner part of the 
disc and/or the primary star. Our preliminary anal- 
ysis indicates a mass ratio (M2/M1) between 0.2-0.4 
for OGLE LMC-SC8-125836 (M08), OGLE LMC- 
SC13-156235 and V393 Sco. 

One of the hallmarks of DPVs is the correla- 
tion between orbital period and superperiod. This 
observational feature could be explained identifying 
the superperiod as the time needed for the disc, to 
reach the 3:1 resonance radius. We conjecture that 
in DPVs primary star does not accumulate gas, be- 
cause it probably rotates at critical velocity. Then 
the transferred matter forms a decretion disc around 
the primary that starts expanding inside its Roche 
lobe. While longer is the orbital period, larger is this 
Roche lobe and longer the superperiod. When the 
disc reaches the 3:1 resonance radius (ry) it starts 
to precess (Whitehurst & King 1991, Lubow 1992; 
WK91 and L92 respectively) and matter starts to 
flow outside the binary system. In our view DPVs 
are semi-detached binaries with a low mass ratio al- 
lowing the formation of a circumprimary disc. 

We make some estimates using the theory of 
decretion discs. In this theory, the disc expands due 
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to viscous stress. The diffusion veiocity at a radius 
r is given by: 

v d (r) » f » a/?VGMTr^ 

(Tutukov & Pavlyuchenkov 2004), where 7 is 
a parameter reiated to the mean free paths of 
turbulent elements and a and (3 are free pa- 
rameters describing a diffusive disc. The time 
needed for a mass element to drift from the star to 
the 3:1 resonance radius (r>) is identified with the 
superperiod and it is given at order of magnitude by: 

Plong ~ T t d = « 77(0.460)^ 

where r\ = \/{afi^JGM\) and we have used the 
approximation r r « 0.46a, where a is the binary 
separation (WK91). 

In the case 7=1, i.e. when the mean free paths 
of turbulent elements in the disc are proportional 
to the distance from the central gravitational object 
(and also proportional to the thickness of the disc, 
Tutukov & Pavlyuchenkov 2004), we get: 

Plong ~ 0.317/a 3 / 2 ~ faffl firli 

this expression gives the right linear dependence 
with the orbital period (P or b)- The proportional- 
ity factor depends on the viscosity properties of the 
disc, parametrized in a and (3 in the Tutukov & 
Pavlyuchenkov 's picture. We need a/3 ~ 1.5 x 10 -3 
to reproduce the order of magnitude of the observed 
DPV correlation. Incidentally, this number is of the 
same order of magnitude that the figure found for 
accretion discs of cataclysmic variables during their 
hot states (Tutukov & Pavlyuchenkov 2004). 
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Fig. 1. Left panel: CMD of 59 LMC DPVs shown as 
filled brown circles, open circles are used for 55 SMC 
DPVs. Additional background (small dots) is a random 
sample of stars in the LMC. We have adopted a dis- 
tance moduli shift -0.5 mag for SMC members. Right 
panel: Graphs showing the P or b-Piong relationship for 
galactic (open squares) and MC DPVs. Note the dif- 
ferences between the populations of the three different 
galaxies (Kolaczkowski & Mennickent, in preparation). 
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Fig. 2. Selected data for the galactic DPV V393 Sco. 
Photometric variability of both periods P or b~ 7.7126 d 
(left up) and Pion 3 = 252.0 d (left down) is shown based 
on ASAS-3 data (Pilecki and Szczygiel 2007). Part of 
one high resolution spectrum showing the Ha line region 
(right, Mennickent et al. in preparation). 



